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Abstract 


We present an alternative method to calculate cross sections for multi-parton scat¬ 
tering processes in the Standard Model at leading order. The helicity amplitudes are 
computed using recursion relations in the number of particles, based on Dyson-Schwinger 
equations whereas the summation over colour and helicity conhgurations is performed by 
Monte Carlo methods. The computational cost of our algorithm grows asymptotically as 
S"", where n is the number of particles involved in the process, as opposed to the n!-growth 
of the Feynman diagram approach. Typical results for the total cross section, differential 
distributions of invariant masses and transverse momenta of partons are presented and 
cross checked by explicit summation over colours. 
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1 Introduction 


The simultaneous production of a large number of energetic partons in high energy col¬ 
lisions of hadrons and leptons, at the TeVatron or, in future, at the LHC and the e+e” 
Linear Collider, gives rise to events with many jets in the hnal state. Quite often, these 
multi-jet events offer an important probe of new physics as for example in the case of 
heavy particle decays in the Standard Model and its extensions, e.g. the Minimal Super- 
symmetric Standard Model. A particularly well known example is the Higgs boson decay 
into four jets through W/Z pairs. Moreover, the multi-jet events provide a signihcant 
background to the discovery channels. The description of these processes is difficult even 
at leading order, because the corresponding amplitudes have to be constructed from a 
very large number of Feynman diagrams, making automation the only solution. As an 
example, in Tab. [H the number of Feynman diagrams relevant for the calculation of 
gg ng is collected. As can be seen, it grows asymptotically factorially with the number 
of particles. 

Table 1: The number of Feynman diagrams contributing to the total amplitude for gg —> 
ng. 


Process 

NpG 

99 2^ 

4 

99 ^9 

25 

99 ^9 

220 

99 % 

2485 

99 Q9 

34300 

99 79 

559405 

i 

00 

10525900 

99 ^9 

224449225 

99 10^ 

5348843500 


Another aspect of dealing with multi-particle amplitudes is the systematic organisation 
of the summation over helicity conhgurations and the SU{Nc) colour algebra. If summa¬ 
tion were performed directly then x 3”^ helicity conhgurations and S”" x x 3”'? 
colour conhgurations would have to be considered, where Ug, n^, nq is the number of 
gluons, quarks and antiquarks respectively while ni stands for the number of fermions 
and massless bosons and n 2 is the number of massive vector bosons. Many of these con¬ 
hgurations do not contribute to the amplitude. However, it is very hard to predict which 
should be kept in advance. The only way to achieve the required efficiency is to use Monte 
Carlo techniques. A further complication is connected to the integration over the multi¬ 
dimensional phase space. In this report we present an approach for efficient tree level 
calculations of matrix elements for multi-parton hnal states which addresses the above 
problems and therefore improves the currently available techniques [1-4]. In this paper. 
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the development of a Monte Carlo summation over colours for the full amplitude and of 
an efficient multi-particle phase space integrator are presented. All algorithms have been 
implemented in a Fortran 95 program that will be subsequently made publicly available. 

The layout of the paper is as follows. In Section [2] the current status of available 
methods is briefly reviewed. Section [3] describes the colour flow decomposition. Section |4] 
presents the recursive relations and algorithms used to build the amplitude. In section [5] 
we give the details of the internal organization of our implementation of the algorithms. 
Here, we also introduce a new approach to the colour structure evaluation. Its compu¬ 
tational complexity is also briefly analyzed. In Section [ 6 l numerical results for the cross 
sections are presented together with distributions of the invariant mass and transverse 
momentum. The hnal section contains our summary and an outlook on future improve¬ 
ments of the algorithm. An Appendix contains a new algorithm for efficient phase space 
point generation. 


2 Dual amplitudes and colour decomposition 

For generality we consider n-gluon scattering 

g{pi, £ 1 , ai) g{p 2 , £ 2 , 02 ) ^ g{p3, £ 3 , 03 )... g{pn, £n, a„) ( 1 ) 

with external momenta {pi}i, helicities {£i}i and colours {aj}” of gluons i = 1 ,... ,n in 
the adjoint representation. The total amplitude can be expressed as a sum of single trace 
terms: 




IeP(2,...,n) 


( 2 ) 


where o'i{2 : n) represent the J-th permutation of the set { 2 ,..., n} and 
represents a trace of generators of the SU{N^) gauge group in the fundamental represen¬ 
tation. For processes involving quarks a similar expression can be derived [5]. One of the 
most interesting aspects of this decomposition is the fact that the Ai{{pi}i, {£i}i) func¬ 
tions (called dual, partial or colour-ordered amplitudes), which contain all the kinematic 
information, are gauge invariant and cyclically symmetric in the momenta and helicities 
of gluons. The colour ordered amplitudes are simpler than the full amplitude because 
they only receive contributions from diagrams with a particular cyclic ordering of the ex¬ 
ternal gluons (planar graphs). For some processes up to six external partons simple and 
compact expressions exist in literature [5-11]. Moreover, for some special helicity combi¬ 
nations, short analytical forms, called the Parke-Taylor helicity amplitudes or Maximally 
Helicity Violating (MHV) amplitudes, are known for general n. They were first obtained 
by the Parke and Taylor [12] and later on proved to be correct in a recursive approach by 
Berends and Giele [13]. Recently, analytical expressions have also been obtained for other 
helicity configurations [14,15], for an arbitrary number of gluons. However, no analytic 
expressions for all helicity configurations are known, and, with the exception of MHV 
amplitudes, the known analytical expressions are usually cumbersome. The simplicity of 
MHV amplitudes suggests that they can be used as the basis of approximation schemes. 
Nevertheless, for large n, the computation of scattering processes is still problematic and 
time consuming. For example to evaluate the full amplitude, the 2”“^ x (n — 1)! conhg- 
urations of colour ordered amplitudes, have to be considered, where 2 ”“^ corresponds to 
the number of helicity conhgurations for massless particles. To obtain the cross section 
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from the n-gluon amplitude one has to square and sum over helicity and colour of the 
external gluons. The squared matrix element can be computed by 

i: = (3) 

where the (n — 1)! x (n — 1)! dimensional colour matrix can be written in the most general 
form as follows: 

C/j = ^ (4) 

l...Vc 

Needless to say that the evaluation of this matrix, is by its own a formidable task in 
the standard approach. 

An important step in the direction of simplihcation of these calculations has already 
been taken by using helicity amplitudes and a better organisation of the Feynman dia¬ 
grams [5,10,11,16-19]. A signihcant simplihcation in these calculations has been made by 
introducing recursive relations [13,20], which express the n-parton currents in terms of all 
currents up to (n — 1) partons. They are based on smaller building blocks which are just 
colour ordered vector and spinor currents dehned for the off mass shell particles. Extend¬ 
ing the recursive approach beyond colour ordered amplitudes, he. to the full amplitude 
in any held theory, is possible. Another approach in this direction based on the so called 
ALPHA algorithm [21,22] or the Dyson-Schwinger recursion equations has been devel¬ 
oped [23-32], where the multi-parton amplitude can be constructed without referring to 
individual Feynman diagrams. In the latter case, apart from the summation over colour in 
the colour how basis, which will be briehy described in the next section, integration over 
a continuous set of colour variables (as well as havour) was introduced. This integration 
technique, however, does not give a straightforward solution for the efhcient merging of 
the parton level calculation with the parton shower evolution. 


3 The colour flow decomposition 

First, let us briehy review the colour approach used in the original version of HELAC [30,31], 
a multipurpose Monte Carlo generator for multi-particle hnal states based on the Dyson- 
Schwinger recursion equations. As was already mentioned, the colour connection or colour 
how representation of the interaction vertices was used in this case. This representation 
was introduced for the hrst time in [33] and later studied in e.g. [31,34]. The advantage 
of this color representation, as compared to the traditional one is that the colour factors 
acquire a much simpler form, which moreover holds for gluon as well as for quark am¬ 
plitudes, leading to a unihed approach for any tree-order process involving any number 
of coloured partons. Additionally, the usual information on colour connections, needed 
by the parton shower Monte Carlo, is automatically available, without any further cal¬ 
culation. In this approach, the gluon held represented as A“ with a = 1,..., — 1 is 

treated as an N^. x traceless matrix in colour space [29]. The new object {A^)ab where 
A,B = 1,..., Aj, can be obtained by multiplying each gluon held by the corresponding 
matrix as follow: 

JV2-1 

{A^)ab = X! ( 5 ) 

a=l 


4 


( 6 ) 


The colour structure of the three gluon vertex is given now by, see Fig. [T]; 

= --{Sad^cf^eb — ^af^cb^ed) 

ai ^ 

where, on the right hand side only products of 5’s appear. This colour structure shows 
how the colour flows in the real physical process, where gluons are represented by colour- 
anticolour states in the colour space, and reflects the fact that the colour remains un¬ 
changed on an uninterrupted colour line. For a four gluon vertex the following expression 




Figure 1: 


Colour flows for the three gluon vertex. 


A D 


B 


G 



C 


F 


H E 

A F 


B 


C 



E 


H 


D C 



A D 



C 


H 


C 


D 


Figure 2 : Colour flows for the four gluon vertex. 
has to be considered: 

E £a\a2X £xasa4,^ai j.a2 J.a3 j.a4 f^\ 

J J ^AB^CD^EF^GH I'J 

ai 

with three permutations of the 010203 indices, which correspond to the six colour flows 
presented in Figl2l Because gluons have different colour states, described for Nc = 3 
within the U{3) group we have an additional unphysical neutral U{1) gluon. This neutral 
gluon does not couple to other gluons as can be easily seen from Eq.([ 6 ]). It couples only 
to quarks and acts as a colourless particle, see second part of Eq.(IH]). The colour structure 
of the quark-antiquark-gluon vertex can be described as follows, see FigJSj 

J2^ABtcD = A^AohcB - -^5 ab5cd)- ( 8 ) 
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Figure 3: Colour flows for the qqg vertex. 

Let us introduce a more compact notation and associate, to each gluon, a label {i, afl 
which refers to the corresponding colour index of previous equations, namely 1 —*■ A, 
(Ji —> -B and so on. The use of labels will be explained latter on. With this notation the 
hrst term of the three gluon vertex is proportional to: 

*^2(T3 • (9) 

For the same graph with inverted arrows a minus sign and interchanged 2 <->■ 2 has to be 
included as well. The momentum part of the vertex, is still the usual one and in 

our notation simply given by: 

g^^Pi - P2? + 9 ^\p2 -PzY + 9^\p, - Pif. (10) 

For the qqg vertex we associate a label (i, 0) for quark and (0, cTj) for antiquark. Finally 
the four gluon vertex is given by a colour factor proportional to: 

^Icra ^3(72 ^2cr4 ^4cri (11) 

with six possible permutations, and a Lorentz part, 

- g^^g^^ - g^^g^^ (12) 

where all three permutations should be included. 

To make use of the colour representation described so far, let us assign, to each external 
gluon, a label (i,cr/(i)), to a quark (i,0) and to antiquark (0,cr/(z))), where i = l...n 
and ai{i),I = 1... n! being a permutation of {1... n}. Since all elementary colour factors 
appearing in the colour decomposition of the vertices are proportional to 6 functions the 
total colour factor can be given by 

«(}?)= E (13) 

/=P(2,...,n) 

^1(7j(1)^2(7j{ 2) • • • ^ncrj{n) j (1^) 

The colour matrix dehned as 

C„ = E»M (15) 

1.J 
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with the summation running over all colours, 1,..., has a very simple representation 
now 

Cjj = (16) 

where 1 < m{ai, aj) < n counts how many common cycles the permutations aj and 
CTj have. The practical implementation of these ideas is straightforward. Given the 
information on the external particles contributing to the process we associate colour labels 
of the form {i,(Ti) depending on their flavour. According to the Feynman rules the higher 
level sub-amplitudes are built up. Summing over all n\ colour connection conhgurations, 
where n is the number of gluons and qq pairs in the process using the colour matrix Cjj 
we get the total squared amplitude. It is worthwhile to note that summing over all colour 
conhgurations is efficient as long as the number of particles is smaller than 0(8). If the 
number of gluons and/or qq pairs is higher than 0(8), since the number of colour hows 
in general grows like (n — 1)!, it also starts to be problematic from the computational 
point of view. For multi-colour processes other approaches have to be considered. The 
natural solution would be to replace the summation over all colour connections by a 
Monte Carlo. However, Monte Carlo summation is not straightforward in the color how 
approach because of the destructive interferences between diherent colour hows that can 
give a negative contribution to the squared matrix element. 

In the limit N,, —>■ cx), only the diagonal terms, I = J, survive, and all 0{N~‘^) 
terms can be safely neglected both in the colour matrix and in the \Ai\^. The interfer¬ 
ences between diherent colour hows vanish in this limit. In the so called Leading Colour 
Approximation (LCA), the squared amplitude for the purely gluonic case is given by 

E mipt}", u*}?. {«.}?) t -1) E E1-4/^ (17) 

a,e e / 

The term — 1 instead of N^, which of course are equivalent in the limit Nc ^ oo, 
has been kept in order to reproduce the exact results for n = 4 and n = 5. In case 
when qq pairs are present the colour factor — 1) is still the same, however 

n = ng + Uqq in this case, where Ug, Uqq is the number of gluons and qq pairs respectively. 
This simplihcation of the colour matrix speeds up the calculation. Moreover, Monte Carlo 
summation over colours can now be performed. 

4 Dyson-Schwinger recursion relations 

Let us now review the recursive relations based on Dyson-Schwinger equations for the cal¬ 
culation of partial amplitudes [32]. These equations give recursively the n—point Green’s 
functions in terms of the 1—, 2—,..., [n — 1)—point functions. They hold all the informa¬ 
tion for the helds and their interactions for any number of external legs and to all orders in 
perturbation theory. We will concentrate here on the gluon, quark and antiquark recursion 
relations, however, in the same way recursive equations for leptons and gauge bosons can 
be obtained. The diagrammatic picture behind the recursive relations is actually quite 
simple. The tree-level recursive equation can be diagrammatically presented as shown in 
FigJU Let Pi,P 2 , ■ ■ ■ ,Pn represent the external momenta involved in the scattering process 
taken to be incoming. In order to write down the recursive relation explicitly, we hrst de- 
hne a set of four vectors [A^(P); (A, B)], which describes any sub-amplitudes from which 
a gluon with momentum P and colour-anticolour assignment A, B can be constructed. 
The momentum P is given as a sum of external particles momenta. Accordingly we dehne 
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a set of four-dimensional spinors ["^(-P); (^,0)] describing any sub-amplitude from which 
a quark with momentum P and colour A can be constructed and by ['^(-P); (0,i?)] a set 
of four-dimensional antispinors for antiquark with anticolour B. The Dyson-Schwinger 



Figure 4: Recursive equation for an off mass shell gluon of momentum P. 


recursion equation for a gluon can be written as follows: 

n 

mpy, (A, B)] = 1:1 A>‘(py-, (A, B).] (18) 

i=l 

+ ^ [ (tg) VP'^\P,py,p 2 )AMMP 2 MpuP 2 ); iA,B) = {C,D), 0 

P=P1+P2 

H [ (/) G""^^^(P,Pl,P2,P3)4li.(pi)4lA(p2)4lp(p3)c^(Pl,P2 +P 3 ); 

P=Pl+P2+P3 

{A,B) = {C,D),^{E,F)2^{G,H)s] 

+ H [ (*^) V'(pi)7'''0(p2)cr(pi,P2); (4l,P) = (0,P>)i (g) (C',0)2] 

P=Pl+P2 

where A, B, C, D, E, E,G,H = 1, 2, 3. The rules for merging colour and anticolour of the 
particles will be explained in the next section. The V^'^^{P,pi,p 2 ) and {P,pi,p 2 ,Ps) 
functions are the three- and four- gluon vertices presented in the previous sections and 
the symbol u(pi,p 2 ) is the sign function which takes into account the Fermi sign when 
two identical fermions are interchanged. The exact form of this function can be found in 
Ref. [30,32]. The sums are over all combinations of pi,P 2 or pi,P 2 ,P 3 that sum up to P. 
The propagator of the gluon is given by: 


.pr _ igpu 
i-i-pu — p2 

For a quark of momentum P we have, Fig IS) 


(19) 



Figure 5: Recursive equation for an off mass shell quark of momentum P. 


















Figure 6: Recursive equation for an off mass shell antiquark of momentum P. 


WP); (A, 0)1 = 5:1 ^(p,)i {A, 0)i] 

2=1 


( 20 ) 


+ H I (*9) P'4'‘(p,)7„^(p2)ff(pi,P2); (A.O) = (B,C)i 

P=P1+P2 


where V is the propagator 


Finally for an antiquark, FiglHl 


V = 


p2- 


(0, ^)] = ^p=Pi V'fe); (0, A)i] 

2=1 


(A 0 ) 2 ! 


( 21 ) 


( 22 ) 


+ Y. I (*9)'/'(P2)2l'‘(pi)7 „'P(t(pi,P 2 ); (O.A) = (B,C)i ® (0, £ 1 ) 2 ] 

P=pi+p2 

where 

P = ^- (23) 

In the same spirit, the recursion equations for all leptons and gauge bosons can be written 
down. With this algorithm we can thus compute the scattering amplitude for any initial 
and hnal states taking into account particle masses as well. In particular any colour 
structure can be assigned to the external legs. Finally this approach has an exponential 
growth of computational time with the number of external particles instead of the factorial 
growth when Feynman graphs are considered. It can be farther optimised in order to 
reduce computational complexity by replacing each four gluon vertex by a three particle 
vertex by introducing the auxiliary field represented by the antisymmetric tensor 
This new field has a quadratic term without derivatives and, therefore, has no independent 
dynamics. The part of the QCD Lagrangian that describes the four gluon vertex 


C = 



Ff. = d,Al 


duAl 


+ 9f 


abc Ab 




(24) 


can be rewritten in terms of the auxiliary held as follows: 


C 


2 / 4(2 ' ^ flU 


(25) 


A single interaction term of the form is left instead of interaction terms rep¬ 

resented by Eq. fIMl) . The recursion for the gluons changes slightly, in fact only for the 
four-gluon vertex part. However, we have an additional equation for the auxiliary held: 


(A, B)] = El (A, B)i] (26) 

2=1 
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P=Pl+V2 


+ Y. I (*9) n; 9l„(Pi)ffA,(P2)ff(pi. P 2 ); (91, B) = (C, £>)i 0 (E, F) 2 ] 


P=P1+P2 


+ H [ (*^) V'(pi) 7'''0 (p 2 )ct(pi,P 2); (94,5) = (0,D)i 0 (C, 0)2] 


^’=Pl+P2 


and 


[ff'“'(F);(9l,S)]= 1 : [(jg) (9'‘V'’-g‘'V'’) A(pi)9l,(P2)<9(pi.P2); (27) 


P=pi+p2 


(24,i?) = (C^,^)i®(9E,9F)2]. 


A few comments are now in order. First, all calculations were performed in the light cone 
representation and all momenta were taken to be incoming. Second, this new held has 
six components. Additionally, as we can see in Eg. (1261) and Eg. ([27]) the colour structure 
of this new vertices, remains the same as in case of the three gluon vertex. The number of 
the types of the sub-amplitudes one has to calculate is in general doubled, however their 
structure is much simpler, which saves computational time while the iteration steps are 
performed. 

After n — 1 steps, where n is the number of particles under consideration, one can get 
the total amplitude. The scattering amplitude can be calculated by any of the following 
relations, depending on the process under consideration. 


where i corresponds to gluon 
where i corresponds to guark line (28) 

'ijj{pi)'ijj{Pi) where i corresponds to antiguark line 


where 


P^ = Y.PP 


so that Pi + Pi = 0. The functions with hat are given by the previous expressions except 
for the propagator term which is removed by the amputation procedure. This is because 
the outgoing momentum Pi must be on shell. The initial conditions are given by 


= eA(!>i).2'= ±1,0 



(29) 


where the explicit form of e'^,ux,vx,ux,vx is given in the Ref. [30]. 

5 Organisation of the calculation 

The recursive algorithm presented in the previous section has the advantage that any 
colour representation can be used in order to assign colour degrees of freedom to the 
external legs in the process under consideration. It can be used either for colour ordered 
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amplitudes or for the full amplitudes as well. The latter case is iu fact the topic of this 
sectiou. Moreover, au alteruative method for takiug iuto accouut the colour structure 
of scatteriug partous, based ou regular colour couhguratiou assigumeuts as compared to 
colour flow oues, will be iutroduced. First, however, the orgauisatiou of the calculatiou 
will be explaiued briefly iu order to better describe the geueral structure. Coutrary to 
the origiual HELAC [30,31] approach, iu this uew versiou the computatioual part cousists 
of oue phase only. This is uot optimal for electroweak processes with moderate uumber 
of exterual particles, but it becomes quite efficieut for processes with mauy particles, 
especially wheu ouly a few species of particles are iuvolved (scalar amplitudes, gluou 
amplitudes iu QCD, etc). The vertices described by the Staudard Model Lagraugiau are 
implemeuted iu fusiou rules that dictate the way the subamplitudes, at each level of the 
recursiou relatiou, will be merged, iu order to produce a higher level subamplitude. Iu 
case wheu a quark is combiued with au autiquark, for example, there are three possible 
subamplitudes which describe three differeut iutermediate states iucorporated iuside the 
fusiou rules, uamely 7 , Z aud g. 

Let us uow preseut the colour mergiug rules which are evaluated iteratively at the 
subamplitude level. Duriug each iteratiou wheu two particles are combiued their corre- 
spoudiug colour assigumeuts are combiued. We have three possibilities to obtaiu a gluou 
{A, B) described by the recursive relatiou Eq. fl26l) . First, it cau be obtaiued wheu a quark 
with colour assigumeut (C, 0) is merged with au autiquark of auticolour (0,T)). Secoud, 
wheu two gluous (C, D) aud {E, F) are combiued, where A, B, C, D,E,F = 1, 2, 3, aud 
dually wheu a gluou aud au auxiliary held H are combiued. Iu the last case, the colour 
structure of the vertex is ideutical to the three-gluou vertex. So we eud up with the 
followiug rules: 


(Ai?)^(C,0)®(0,Z}), (30) 

{A,B)^{C,D)^{E,F). (31) 

Iu the hrst case the mergiug occurs accordiug to the rule preseuted iu Eq.([H]). aud the 
followiug gluou cau be produced: 

(A5) = (c,0)® (0,D) = (C,/!), iiC^D. (32) 


However, the situatiou is more complex wheu a quark aud au autiquark have the same 
colour aud auticolour: 

(H, B) = {C, 0) 0 (0, D) = (1,1)^, © (2, 2)^, © (3, 3 )^ 3 , A C = D. (33) 

with = 0. We thus have three possibilities with differeut weights, which for iustauce 
are Wi = W 2 = —1/3, = +2/3 iu case C = D = 3. 

Iu the secoud case, wheu gluous are combiued accordiug to Eq.([ 6 ]) we have the followiug 
optious: 


{A,B) = {C,D)(^{E,F) = {E,D), if C = F aud E 7 ^ D (34) 

{A,B) = {C,D)(Z{E,F) = {C,F), if D = E and C ^ F. (35) 

Moreover, we have au additioual possibility wheu colours aud auticolours of gluous are 
the same. Oue gets, iu this case, a gluou iu two colour states with differeut weight: 

(H,H) = (O,T)) 0 (E,F) = (0,0)+i©(T),T))_i, if C = F E = D aud C ^ D. (36) 
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Finally for C = D = E = F the result vanishes identically. 

To obtain the quark, (A, 0), described by the recursive relation Eq. fl20|) we have to 
combine a gluon with another quark one more time according to the rule presented in 
Eq.®: 

(AO) = (5,C'){8)(D,0). (37) 

We have two possibilities for colour assignment: 

(AO) = (5,C')®(D,0) = (5,0), ifC = 5, (38) 

(71,0) = (5,(7) ® (5,0) = (5,0), if 5 = 5. (39) 

For the antiquark the situation is the same, so we will not elaborate it here. 

For the qq interactions with 7 , or Higgs boson the situation is very simple and 

we have only one possibility: 

(71,0) 0 (0,5) = (0,0) if 71 = 5. (40) 


The sum over colour can be performed in this way by considering all possible colour- 
anticolour conhgurations according to the above rules. However, the procedure can be 
facilitated by Monte Carlo methods where a particular colour-anticolour conhguration 
is randomly selected. As we will see an important gain in computational efficiency is 
achieved within this framework. 

The necessary condition which must be fulhlled, while the particular colour assignment 
for the external coloured particles is chosen, is that the number of colour and anticolour 
of each type, is the same. Otherwise the particles can not be connected by colour flow 
lines and the amplitude is identically zero. In the Monte Carlo over colours method one 
has to multiply the squared matrix element by a coefficient that counts the number of 
non zero colour conhgurations. 

The number of non zero colour conhgurations, according to the above mentioned 
necessary condition is given by 


Ncc 


riq riq—A 

E E 


A=0 B=0 


^0 


A + B + C) 


(41) 


where Ug is the total number of colours and A, B,C are the numbers of the colour type 1 , 
2 and 3 respectively. As we already stated, the condition given by Eq. fl4T]) is necessary but 
not sufficient. Among this set there are still conhgurations which do not give contributions 
to the total amplitude. 

In Tab. [2l and Tab. O diherent numbers of colour conhgurations in the process 
under consideration are presented. In the hrst column the total number of all colour 
conhgurations is listed, where and Ug, Ug are the number of quarks and 

antiquarks respectively and gluons are treated as qq pairs. The second column represents 
results for the number of non vanishing colour conhgurations N^c calculated using Eq. flTT|) . 
In the next column the ratio is presented. The number of colour conhgurations 

(in percentage) inside the N^c set which hnally gives rise to non zero amplitudes is shown 
in the last column. Those numbers are evaluated by Monte Carlo. Note that while 
the number of external particles is increased the corresponding number of vanishing colour 
conhgurations in the third column is decreased, as we can see in Tab. [2l 

As far as the summation over the helicity conhgurations is concerned there are two 
possibilities, either the explicit summation over all helicity conhgurations or a Monte 
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Table 2: The number of colour configurations for the processes with gluons only. 
corresponds to all possible colour configurations, while N^c corresponds to the number 
of colour configurations calculated using formula Eg. In the third column the ratio 
Ncc/Nqc^ is presented. In the last column the number of non vanishing colour configura¬ 
tions evaluated by MC (in percentage) inside the Ncc is shown. 


Process 

atALL 

Ncc 

Nco/NS" 

N"cc (%) 

99 2^ 

6561 

639 

0.0974 

59.1 

99 ^9 

59049 

4653 

0.0788 

68.4 

99 ^9 

531441 

35169 

0.0662 

77.4 

99 I>9 

4782969 

272835 

0.0570 

85.0 

99 Q9 

43046721 

2157759 

0.0501 

90.4 

99 19 

387420489 

17319837 

0.0447 

94.0 

i 

00 

3486784401 

140668065 

0.0403 

96.4 


Carlo approach. In the latter case, for example for gluon, it is achieved by introducing 
the polarisation vector 

£j(p) = e‘«£^(p) + e-‘«£'l(p), (42) 

where 0 G (0, 27r). By integrating over (j) we can obtain the sum over helicities 

1 r^TT _^ 

— dfi e^(p)(e^(p))* = Y) e^(p)(e’)(p)y. 

A=i 


To determine the cost of computation of the n—point amplitude using the algorithm 
based on Dyson-Schwinger equations, one has to count the number of operations [35]. 
( n\ 

There are I ^ 1 momenta at each level. The total amount of sub-amplitudes corre¬ 
sponding to those momenta is simply given by: 

E ( ) = 2” - 2 (43) 

where n is the number of particles involved in the calculation. Moreover one has to count 
how many ways exist to split a number of level k to two numbers of levels ki and k 2 . 
The last step is to sum over all levels. The total number of operations that should be 
performed in the case when only three-point vertices exist is then 

E ( r ) E ( = E ( r ) {2'' - 2} = 3” - 3.2” + 3. (44) 

In the limit n —>■ cx) the number of operations grows like 3” instead of the n\ growth in the 
Feynman graph approach. When the Monte Carlo over colour structures is performed and 
only one particular colour conhguration is randomly chosen the computational cost of this 
algorithm is given exactly by this expression. Otherwise, this formula must be multiplied 
by the number of non zero colour conhgurations for the process under consideration. 
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Table 3: The number of colour configurations for the processes with gluons and one or two 
qq pairs. corresponds to all possible colour configurations, while N^c corresponds to 

the number of colour configurations calculated using formula Eg. (f^^- In the third column 
the ratio is presented. In the last column the number of non vanishing colour 

configurations evaluated by MG' (in percentage) inside the Ncc is shown. 


Process 


N 


cc 


N, 


Ncc/N; 


ALL 

CC 


NL (%) 


99 

99 

99 

99 

99 

99 

99 


99 

99 

99 

99 

99 


uu 

729 

guu 

6561 

2guu 

59049 

3guu 

531441 

Aguu 

4782969 

5guu 

43046721 

6guu 

387420489 


93 

0.1276 

639 

0.0974 

4653 

0.0788 

35169 

0.0662 

272835 

0.0570 

2157759 

0.0501 

17319837 

0.0447 


93.5 

91.6 

92.6 

94.6 
96.4 
97.8 

98.6 


cccc 

gcccc 

2gcccc 

Sgcccc 

igcccc 


6561 

59049 

531441 

4782969 

43046721 


639 

4653 

35169 

272835 

2157759 


0.0974 

0.0788 

0.0662 

0.0570 

0.0501 


99.1 

98.8 

99.0 

99.3 

99.6 


6 Numerical Results 

In this section numerical results for multi-parton production at the LHC are presented. 
As we will show the Monte Carlo summation over colours, speeds up the calculation 
substantially, compared to the one based on explicit summation, especially for processes 
with a large number of colors, i.e. gg —> ng with n > 6. 

The centre of mass energy was chosen to be ^/s = 14 TeV. In order to remain far from 
collinear and soft singularities and to simulate as much as possible the experimentally 
relevant phase-space regions, we have chosen the following cuts: 


PTi > 60 GeV, \yi\ < 2.5, ARij >1.0 


(46) 


for each pair of outgoing partons i and j. Here pxi and yi are the transverse momentum 
and rapidity of a parton respectively defined as: 


Pt, = 




( Ej -I- pzj \ 

V-Bi -P.-J ' 


(46) 


In practice for massless quarks the rapidity is often replaced by the pseudorapidity variable 
rj = — lntan(0/2), where 6 is the angle from the beam direction measured directly in the 
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Table 4: Results for the total cross section for processes with gluons and quarks with up 
to two qq pairs. (Jexact corresponds to summation over all possible colour configurations, 
while (Jmc corresponds to Monte Carlo summation. 


Process 

f^EXACT A E (nb) 

£(%) 

(Jmc ± £ (nb) 

e(%) 

99 2^ 

(0.46572 ± 0.00258) X10^ 

0.5 

(0.46849 ± 0.00308) X10^ 

0.6 

99 ^9 

(0.15040 ± 0.00159) X10^ 

1.0 

(0.15127 ± O.OOllojxlO^ 

0.7 

99 % 

(0.11873 ± 0.00224) xlO^ 

1.9 

(0.12116 ± 0.00134)x102 

1.1 

99 % 

(0.10082 ± 0.00198)xl0^ 

1.9 

(0.09719 ± 0.00142)xl0i 

1.5 

99 Q9 

(0.74717 ± 0.01490) X10-1 

2.0 

(0.76652 ± 0.01862) X10-1 

2.4 

99 uu 

(0.36435 ± 0.00199) xl02 

0.5 

(0.36619 ± 0.00132)x102 

0.4 

99 9UU 

(0.35768 ± 0.00459) xlOi 

1.3 

(0.35466 ± 0.00291) xlOi 

0.8 

99 ‘^9uu 

(0.49721 ± 0.00758) X10° 

1.5 

(0.50053 ± 0.00725) X10° 

1.4 

99 ^9uu 

(0.50598 ± 0.01441)xl0-i 

2.8 

(0.52908 ± 0.01264) X10-1 

2.4 

99 4:guu 

(0.51549 ± 0.02017)x10-2 

3.9 

(0.51581 ± 0.01245)x10-2 

2.4 

gg —>• cccc 

(0.25190 ± 0.00528) X10-2 

2.1 

(0.24903 ± 0.00373) X10-2 

1.5 

99 9CCCC 

(0.60196 ± 0.01908)xl0-° 

3.2 

(0.58817 ± 0.00926)xl0-° 

1.6 

gg — 2gcccc 

(0.95682 ± 0.03441) X10-^ 

3.6 

(0.92212 ± 0.02485) X10-^ 

2.7 


detector. The last variable is ARij which is the radius of the cone of the parton dehned 
as 

AR,, = + (,,i - (47) 

with azimuthal angle = <hj — G (0 , tt) 




= arccos 


PxiPa 


Py^Pyj 


PTiPTj 


(48) 


Quarks are treated as massless. All results are obtained with a hxed strong coupling con¬ 
stant (q!s=0.13). For the parton structure functions, we used CTEQ6L1 PDF’s parametri- 
sation [36,37]. For the phase space generation we used the algorithm described in the 
Appendix, whereas in several cases results were cross checked with 


• PHEGAS [38], which automatically constructs mappings of all possible peaking struc¬ 
tures of a given scattering process and uses self-adaptive procedures like multi¬ 
channel optimisation [39]. PHEGAS exhibits a high efficiency especially in the case 
of n(-l- 7 )—fermion production in e^e“ collisions as long as the number of external 
particles is smaller than 0(8). 
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Table 5: Comparison of the computational time for the squared matrix element calculation 
for processes with gluons and/or quarks with up to two qq pairs, Iexact denotes the time for 
the processes calculated with summation over all possible colour flows (for colour ordered 
amplitudes), while Imc is obtained with Monte Carlo summation over colours (for full 
amplitudes). In the last column their ratio is presented. Time values are given in some 
arbitrary units. 


Process 

fCF 

^EXACT 

Imc 

tEXACT/tMC 

99 2^ 

0.315x10° 

0.554x10° 

0.57 

99 ^9 

0.329x101 

0.143x101 

2.30 

99 ^9 

0.383x102 

0.372x101 

10.29 

99 5g 

0.517x10° 

0.105x102 

49.24 

99 Q9 

0.987x10^ 

0.362x102 

272.65 

99 uu 

0.260x10° 

0.466x10° 

0.56 

99 9UU 

0.196x101 

0.123x10° 

1.59 

99 ‘^9uu 

0.166x102 

0.348x101 

4.77 

99 ^9uu 

0.171x10° 

0.129x102 

13.25 

99 4:guu 

0.197x10^ 

0.307x102 

64.17 

gg —cccc 

0.697x101 

0.605x101 

1.15 

99 9CCCC 

0.568x102 

0.217x102 

2.62 

gg — 2g cccc 

0.619x10° 

0.401x102 

15.44 


• HAAG [40], which efficiently maps the so called “antenna momentum structures” 
typically occurring in the QCD amplitudes. HAAG becomes less efficient for processes 
with higher number of particles, more than 0(8), as is the case for any multi-channel 
generator, where the computational complexity problem arises when the number of 
channels increases. 

• RAMBO [41], a flat phase-space generator RAMBO generates the momenta distributed 
uniformly in phase space so that a large number of events is needed to integrate 
the integrands to acceptable precision, which results to a rather low computational 
efficiency. 

Total rates and various kinematical distributions are examined. The algorithm de¬ 
scribed in the previous sections as well as in Ref. [42,43] has been used to compute total 
cross sections for many parton production processes. We give the result with summation 
over all possible colour configurations, called ’exact’, cTexact, as well as the result obtained 
with Monte Carlo summation over colours cTmc- 
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Table 6 : Results for the total cross section for processes with gluons, Z, and quarks 
with up to qq pairs for the higher number of external partons. Cmc corresponds to Monte 
Carlo summation over colours. 


Process 

(Jmc ± £ (nb) 

£(%) 

99 7g 

(0.53185 ± 0.01149)x10-2 

2.1 

i 

00 

(0.33330 ± 0.00804) X10-3 

2.4 

99 ^9 

(0.13875 ± 0.00430) X10-^ 

3.1 

99 5guu 

(0.38044 ± 0.01096) X10-3 

2.8 

gg —> Sgcccc 

(0.95109 ± 0.02456) X10-3 

2.6 

gg —>• Agcccc 

(0.81400 ± 0.02583) xlO-® 

3.2 

gg Zuugg 

(0.18948 ± 0.00344) X10-3 

1.8 

gg ^ W+udgg 

(0.62704 ± 0.01458) X10-3 

2.3 

99 ZZuugg 

(0.16217 ± 0.00420)xl0-® 

2.6 

gg —> W~^W~uugg 

(0.27526 ± 0.00752) X10-3 

2.7 

dd — Zuugg 

(0.38811 ± 0.00569)x10-3 

1.5 

dd —> W~^csgg 

(0.18765 ± 0.00453) X10-3 

2.4 

dd —> ZZgggg 

(0.99763 ± 0.02976) X10-^ 

2.9 

dd W+W-gggg 

(0.52355 ± 0.01509) xlO-® 

2.9 


In case of the explicit summation both colour flow and colour conhguration decompo¬ 
sition have been used to cross check results. As far as helicity summation is concerned, a 
Monte Carlo over helicities is applied. 

The results presented for the total cross sections, have been obtained for 10® Monte 
Carlo points passing the selection cuts given by Eg. dTSj) . In the Tab.lH the results for the 
total cross section for processes with gluons and/or quarks with up to two qq pairs are 
presented. All cross sections are in agreement within the error. For the same number of 
accepted events the results with the Monte Carlo summation over colours can be obtained 
much faster and the error is at the same level compared to the ’exact’ results for the same 
processes. 

In the Tab. [5l a comparison of the computational time for squared matrix element 
calculations for processes with gluons and/or quarks with up to two qq pairs is presented. 
In each case Iexact means time obtained for the processes calculated with the summation 
over all possible colour flows (for colour ordered amplitudes), while Imc is obtained with 
Monte Carlo summation over colours (for full amplitudes). 
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In the next Tab. El the results for the total cross section for processes with gluon, 
Z, and quarks with up to qq pairs for a larger number of external partons are also 
presented. 
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Figure 7: Distribution in z = |A^oneP /where |A^oneP is the square matrix 
element for one particular colour configuration normalised to the sum of all possible. The 
left-hand side plot corresponds to the gg —> 2g process, the right-hand side one to the 
gg —> 3g. Solid line crosses denote summation over all colour configurations whereas 
dashed, the Monte Carlo summation. 





Figure 8: Invariant mass distribution of 2 gluons in the multigluon 2g —*■ ng process. 


In order to demonstrate that the Monte Carlo summation over colours does give the 
same information on the colour connection structure of the process, we examine in Fig. [7| 
the distribution of the following variable. 


l-MoneP 

U=i 


(49) 
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Figure 9: Transverse momentum distribution of a gluon in the multigluon 2g —*■ ng process. 



(CeV) 


Figure 10: Transverse momentum distribution of the hardest gluon in the multigluon 
2g —>■ ng process. 

where jAloneP is the square matrix element for one particular colour connection or colour 
flow conhguration, normalised to the sum of all possible ones. In case of g{l)g{2) 
g{3)g{4:) process, the following colour flow (1—s>3—*>4—>2—>-1) has been plotted. 
This chain of numbers shows how gluons are colour connected with each other. For the 
g{l)g{2) — g{3)g{A)g{5) process in the Fig. [71 the colour flow (1—>5—>3—*>2—>4—*>1) 
is used. 

The 2 ;—variable distributions can be used in order to extract colour connection infor¬ 
mation needed by a parton shower calculation on the event by event basis. The agreement 
between the ’exact’ and MC distributions implies that we will get the same information 
on the colour structure of the amplitude and that the merging of parton level calculation 
with the parton shower evolution can be safely performed in order to achieve a complete 
description of the fully hadronised final states observed in real experiments. 

Finally, the distribution of the invariant mass of 2-partons, the transverse momentum 
distribution of one parton as well as of the hardest and the softest parton for gg ^ ng, 
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Figure 11: Transverse momentum distribution of the softest gluon in the multigluon 2g —> 
ng process. 

u = 4,..., 8 process is shown in Fig. [HI— Fig. [11] Distributions of the Monte Carlo over 
colours clearly demonstrate that this approach performs very well not only at the level of 
total rates but also at the level of differential distributions. 


Summary and Outlook 

In this work an efficient way for the computation of tree level amplitudes for multi-parton 
processes in the Standard Model was presented. The algorithm is based on the Dyson- 
Schwinger recursive equations. We discussed how the summation over colour conhgura- 
tions can be turned into a Monte Carlo summation, which proved to be more efficient, 
especially for a large number of coloured partons. Additionally, a set of typical results 
for total cross sections and differential distributions has been given. Moreover, a new al¬ 
gorithm for phase-space generation has been presented and used. The complete package 
can be used to generate, efficiently and reliably, any process with any number of external 
legs, for n < 12, in the Standard Model. 

Our future interest includes a systematic study of fully hadronic hnal states in pp and 
pp collisions which requires the merging of the parton level calculations with parton shower 
and hadronization algorithms, e.g interfacing our package with codes like PYTHIA [44] or 
HERWIG [45]. The development of this kind of multipurpose Monte Carlo generators will 
certainly be of great interest in the study of TeVatron, LHC and e’''e“ Linear Collider 
data. 
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Appendix 

In case of scattering of two hadrons it is useful to describe the final state in terms of 
transverse momentum azimuthal angle 0 and rapidity y variables. These variables 
transform simply under longitudinal boosts which is useful in case of the parton-parton 
scattering where the centre of mass system is boosted with respect to that of the two 
incoming hadrons. In terms of 0 and y the four-momentum of a massless particle can 
be written as 


p^ = {E, p^, p^) = {pt cosh p, pt cos 0, pr sin 0, px sinh y) (50) 

It is much more natural to express the phase space volume for a system of n particles 

„ n n 

Vn= 6^{P - X^Pi) n - "l-)0(p-) (51) 


2 = 1 2=1 


using Pt, 0 and y variables 

« n n 

14 = / S‘(P - YlPi)WPT,iPTiiyiiA- 


(52) 


2=1 2=1 


To derive the above expression, which lead us to the Monte Carlo algorithm, we follow 
the method presented in Ref. [41]. We start by defining the phase-space-like object 


Rn = 


[Q dkTi Pikq 


\i=2 


n27T ( P' 

\\{d(t)i 
^ \i=2 


^+oo / ^ 

n 

\i=2 


(53) 


describing a system of n four momenta that are not constrained by momentum conser¬ 
vation but occur with some weight functions P{kTf), n(pj) which keeps the total volume 
hnite. In the next step we have to relate the new variables to the physical ones pjv, p* 
and 0 *: 

n-\-oo 


14 = 


( I n *T, I ^ 

^ n ^ 

n dPT^ - xkr^) 


r2TT 


Y[d(j)i 




\j=2 

C+oo 


Ki=2 


\{dyi 5{yi + yi_i - yi) 


\i=l 


\i=2 


poo pZTT 

/ dkx., / d 0 i S{xY] kTi cos cfi) kxi sin 00 Ji 

•^0 -^0 *=2 *=2 

poo p-\-oo ^ ^ 

/ dx dyi 5(x V kxi cosh yi - E) (5(a; V kxi sinh p^ - L) J 2 . 

JO J-oo 


(54) 


with the Jacobians Ji and J/ 2 : 

d{x E kxi cos 0i, X E kxi sin 00 


Ji = 
J2 = 


0(/CTi,00 


= x^k 


Ti 


d{x E kTi cosh yi - E,xJ2 kxt sinh yi - L) 


<9(x,p0 


^2 


X 


(55) 

(56) 
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E and L represents the energy and longitudinal parts of the initial two particles. We 
proceed with integration where the different arguments of the various 5 functions were 
manipulated in order to perform the integral 


/o 


2 = 1 
+ 00 ^ 


-xkr,) 

(57) 

Pi-i - Pi)- 

(58) 


' —00 


2=2 


We are left with 


Vo= f S^iP - (f[PTidpTidyid(pi] (f[P - L^)‘^dx. (59) 

i=l \i=l / \i=2 \ X J PTi ) X 

We are free to choose the distribution functions so that the total volume is kept hnite. 
The criterion used is to minimize the variance, by taking into account the anticipated 
form of the multi-parton matrix elements, so we introduce 


1 


—X 


P{x) = -exp {— , n(?/) = 


tanh( 2 r 7 + y) + tanh( 2 r 7 — y) 


8p 


where a > 0 and perform the integration over dx 


^ n 1 

n-exp 


\i=2 


PTi\ 1 \ 1 


—n +1 


xa j pTi 7 X' 


Ki=2 PXi / \i=2 


We hnally arrive at the formula 

P n / n \ 

Vo = 5^{P - Y.Pi) \ X{PTidpTidyid(l)i 


X 


2 = 1 \2 = 1 


—n+l 


IS'’" 


r(n - 1 ) n n(s,) (-B" - i") 2 ' 


\i=2 


On the other hand if we applied Eq. fl60|) to the formula Eq. flS^ we can find 


dkTi - exp 


-k 


T- 

X t 


= 1 


\i=2 


r+oo / \ 

(l[dmy^)] = 

\i=2 / 


10 


2Tr / ^ 

E 

i=3 


n #. =(2 


TT 


\n—1 


and 


w = (27r ■ 8pr-\ 


The weight of the event is given by 


W = 


{2tt ■ 8r/)"-i 


(60) 


= n — Y.Pt, r(n - 1 ), ( 61 ) 


(62) 


(63) 

(64) 

(65) 

( 66 ) 

(67) 
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where 


s„ = I dpT.dyd<Pi{ujA (^Pt) r(n - 1) ("h nfe) Ve" - L") 2" 

In the next step we translate this description into a Monte Carlo procedure and generate 
independently n variables /ct;, 4>i and yi = yt — yi-i and assuming that 0 i = 0 as well as 
yi = 0. Using the symbol pi to denote a random number uniformly distributed in (0,1) 
we do this as follows: 


kxi = -a log Pi, i = 2,...,n ( 68 ) 

(pi = 271 Pi, i = 2,...,n (69) 

where a is a free parameter. For yi variable we proceed in few steps starting with 

Fi = exp(4p(2pi - 1)) i = 2,...,n (70) 


where p is a free parameter and 


, (Fi + 1 ) sinh 2 p 
cosh Pi =-—-, 


. , _ (Fi —l)cosh2p 
smhpi =-—- 


(71) 


where 

Zi = ^2F,cosh4p-(l + F;2)_ (72) 

To complete the description of the algorithm we have to find an expression for /c^i, (pi and 
Pi variables. We start by defining the transversal part of the 2,... ,n system as follows 


n n 

^ = ^^Ti COS (pi, ^ sin (pi. 

i=2 i=2 


(73) 


/^From these equations we have the following relations for and (pi to be able to describe 
the total n particle system 


cos (pi 


X 


sin 01 = 


Y 

kxi 


where 

A)Ti = + Y^. 

The total energy and total longitudinal part of the system are represented by 


(74) 

(76) 


S = hr, + kr-i cosh(y2) + kr, cosh(!/s + 1 / 2 ) + ■ ■ ■ 
C ^ sinh(y2) + sinh(y3 + 112) + ■ ■ ■ 

SO Pi is given by 

SE-CL 

v/£2 _ 

-CE + £L 


cosh Pi = 
sinh Pi = 


(76) 

(77) 


(78) 

(79) 


£2 _ 1,2 














Finally to get the final four momenta pf the following transformations are used: 


Ei = xkxi cosh Pi 

(81) 

Pi = xkTi cos (pi 

(82) 

pf = xkr^ sin (pi 

(83) 

Pi = xkTi sinh pi. 

(84) 


This completes the description of the algorithm we have to supplemented it with the 
prescription for the weight of a generated event which is given by Eg . 0671) . 
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